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Abstract. We describe a six-parameter family of the minimum- uncertainty squeezed states for 
the harmonic oscillator in nonrelativistic quantum mechanics. They are derived by the action of 
corresponding maximal kincmatical invariance group on the standard ground state solution. We show 
that the product of the variances attains the required minimum value 1/4 only at the instances that 
one variance is a minimum and the other is a maximum, when the squeezing of one of the variances 
occurs. The generalized coherent states are explicitly constructed and their Wigner function is 
studied. The overlap coefficients between the squeezed, or generalized harmonic, and the Fock 
states are explicitly evaluated in terms of hypergeometric functions. The corresponding photons 
statistics are discussed and an application to quantum optics and cavity quantum electrodynamics is 
mentioned. Explicit solutions of the Heisenberg equations for radiation field operators with squeezing 
are found. 



l. An Introduction 

From the very beginning, nonclassical states of the linear Planck oscillator, in particular the 
coherent and squeezed states, have been a subject of considerable interest in quantum physics (see 
[30] . [38] . [51], [77], [79], |130j . |131] and the references therein). They occur naturally on an atomic 
scale [12], [75] and, possibly, can be observed among vibrational modes of crystals and molecules 
[10], [IE]- A single monochromatic mode of light also represents a harmonic oscillator system for 
which nonclassical states can be generated very efficiently by using the interaction of laser light 
with nonlinear optical media [15], [91], [99], [100], [102], [fl9], [128], [152]. Generation of squeezed 
light with a single atom has been experimentally demonstrated |116j . On a macroscopic scale, the 
squeezed states are utilized for detection of gravitational waves [66] below the so-called vacuum 
noise level and without violation of the uncertainty relation [T], [12], |118j . |148j . 

The past decades progress in generation of pure quantum states of motion of trapped particles 
provides not only a clear illustration of basic principles of quantum mechanics, but it also manifests 
the ultimate control of particle motion. These states are of interest from the standpoint of quantum 
measurement concepts and facilitate other applications including quantum computation (see [TT] . 
PH, [20], [22], [52], [59], [61], [75], [90], PS], [W], [HD], [117], [123], [126] and the references 
therein) . 
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It is well known that the harmonic quantum states can be analyzed through the dynamics of a 
single, two-level atom which radiatively couples to the single mode radiation field in the Jaynes- 
Cummings(-Paul) model [19], [23], [73], [H], p9], [131], [ED] extensively studied in the cavity 
QED [UJ, [52], |12U] . |121j . Creation and detection of thermal, Fock, coherent, and squeezed states 
of motion of a sine le 9 Be + ion confined in a rf Paul trap was reported in |106] , where the state of 
atomic motion had been observed through the evolution of the atom's internal levels (e.g., collapse 
and revival) under the influence of a Jaynes-Cummings interaction realized with the application of 
external (classical) fields. The distribution over the Fock states is deduced from an analysis of Rabi 
oscillations. 

Moreover, Fock, coherent, and squeezed states of motion of a harmonically bound cold cesium 
atoms were experimentally observed in a ID optical lattice [TTj . jllOj . This method gives a direct 
access to the momentum distribution through the square of the modulus of the wave function in 
velocity space (see also [IH],[ig,[5I],[21],[22],[Sl],[Zl],[Z5],[SD], [III], [HH] and the references 
therein regarding cold trapped ions and their nonclassical states; progress in atomic physics and 
quantum optics using superconducting circuits is reviewed in [52], |157j ). 

Recent reports on observations of the dynamical Casimir effect [53], |156j strengthen the interest 
to the nonclassical states of generalized harmonic oscillators [30], [31], [35], [36], [37], [58], [103] , 
|104j . |112j . |144] and [151] . The amplification of quantum fluctuations by modulating parameters 
of an oscillator is closely related to the process of particle production in quantum fields [31], [71] . 
[104J, and [112J. Other dynamical amplification mechanisms include the Unruh effect [147] and 
Hawking radiation [TP] , [52] , [65] . 

The purpose of this paper is to construct the minimum- uncertainty squeezed states for quantum 
harmonic oscillators, which are important in these applications, in the most simple closed form. 
Our approach reveals the "hidden" quantum numbers/integrals of motion of the squeezed states in 
terms of solution of certain Ermakov-type system [91], [95]. The corresponding generalizations of 
Fock states, which were originally found in [105] and recently rediscovered in [95], are discussed in 
detail. As a result, the probability amplitudes of these nonclassical states of motion are explicitly 
evaluated in terms of hypergeometric functions. Their experimental observations in cavity QED and 
quantum optics are briefly reviewed. Moreover, the radiation field operators of squeezed photons, 
which can be created from the QED vacuum, are introduced by second quantization with the aid 
of hidden symmetry of harmonic oscillator problem in the Heisenberg picture. 

In summary, experimental recognitions of the nonclassical harmonic states of motion have been 
achieved through reconstruction of the Wigner function in optical quantum-state tomography [15j . 
[102] . from a Fourier analysis of Rabi oscillations of a trapped atom [106] , and/or by a direct 
observation of the square of the modulus of the wave function for a large sample of cold cesium 
atoms in a ID optical lattice [TTJ, [110] . Our theoretical consideration complements all of these 
advanced experimental techniques by identifying the state quantum numbers from first principles. 
This approach may provide a guidance for engineering more advanced nonclassical states. 

The paper is organized as follows. In sections 2 and 3, we describe the minimum-uncertainty 
squeezed states for the linear harmonic oscillator in the coordinate representation. The generalized 
coherent, or TCS states, are constructed in section 4. In sections 4 and 5, the Wigner and Moyal 
functions of the squeezed states are evaluated directly from the corresponding wave functions and 
their classical time evolution is verified with the help of a computer algebra system. The eigenfunc- 
tion expansions of the squeezed (or generalized harmonic) states in terms of the standard Fock ones 



MINIMUM-UNCERTAINTY SQUEEZED STATES 



3 



are derived in section 6 (see also [34] , [39] and the references therein for important special cases). 
Some experiments on engineering of nonclassical states of motion are analyzed in section 7. Here, 
the experimentally observed probability distributions are derived from our explicit expression for 
the probability amplitudes obtained in the previous section. In section 8, we revisit the radiation 
field quantization in a perfect cavity, which is important for applications to quantum optics. Non- 
standard solutions of the Heisenberg equations of motion for the electromagnetic field operators, 
that naturally describe squeezing in the Heisenberg picture, are found. The variance of the number 
operator, which together with the eigenfunction expansion allows one to compare our results with 
experimentally observed squeezed photon statistics [15], [128] . is evaluated from first principles in 
section 9. A brief summary is provided in the last section. A convenient complex parametrization 
of the Schrodinger group can be found in appendix. 



2. The Minimum-Uncertainty Squeezed States 

The Heisenberg Uncertainty Principle is one of the fundamental laws of nature and the coherent 
states that minimize this uncertainty relation are well known. But, equally important in recent 
developments, minimum- uncertainty squeezed states are not so familiar outside a relatively narrow 
group of experts. Here, for the benefits of the reader, we construct these states as explicitly as 
possible and elaborate on some of their remarkable features. 

The time-dependent Schrodinger equation for the simple harmonic oscillator in one dimension, 

2iiP t + iP xx -x 2 iP = 0, (2.1 

has the following square integrable solution (Gaussian wave packet )Q 

i(a(t)x 2 +8(t)x+K(t)+~f(t)) 

^ (*, t) = G - — — (2 . 2 

y/fl (t) y/ir 

where 



li (t) — /i y/3 sin t + (2a sin t + cos t) , (2.3 
_ a cos2t + sin2t (gg + 4qg-l)/4 
(3 sm t + (2«o sin t + cos t) 

P(t)= , & =, (2.5 

/3q sin 2 1 + (2a sin t + cos t) 

i \ 1 0a tant , 

7 it) = 7o - o arctan TT^> — ; — 7' 2 - 6 
2 1 + 2a taut 



^ ^ _ S (2« sin t + cos t) + e (3l sin t ^ 
(3q sin 2 1 + (2«o sin t + cos t) 2 



e (2a sin t + cost) - f3 5 sin t 
e (t) - ; , [Z.x 

/3q sin 2 1 + (2ao sin t + cos t) 



1 This solution is invariant under the time reversal t — > —t, ip — > tp* with ao — > — cxq, 7 — > — 7oj <^o — ^ — 5o, and 
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K (t) = Kq + sin 2 1 



+ -sin 2t 
4 



/3q sin 2 1 + (2a sin t + cos i) 2 

2 o2 _ s.2 

E oPo °o 



(2.9) 



/3q sin 2 1 + (2a sin t + cos £) 2 

(/i > 0, cto, j5 7^ 0, 7 , Sq, Eq, k are real initial data of the corresponding Ermakov-type system; 
a complex form of equations (I2.3D — ( I2T9T) is provided in Appendix A and the invariants are given by 
f l6.lip -f l6.12p ). This quantum state is the special case n = of a 'nonclassical' oscillator solution 
found in |105j . which has also been recently derived in a unified approach to generalized harmonic 



oscillators (see, for example, [25] 



87] . [95] and the references therein). These solutions are 



verified by a direct substitution with the aid of Mathematica computer algebra system [80] (see also 



SB), 



and [97J. (The simplest special case /x = ±/3 Q = 1 and a = 7 = 5 = e = k = 



reproduces the ground oscillator state [51 



[86] . [108] ; see also the original Schrodinger papers 
[130] . [131] . For the coherent states a = and (3 = ±1, see [131j and a more general case when 
a = is discussed in [57] . [HE]- More details on the derivation of these formulas can be found in 
Refs. [M], [95], and [105] . An analog of Berry's phase is evaluated in Refs. [145] . [146] .) 



The "dynamic harmonic oscillator ground state" fl2.2p - fl2.9l) is the eigenfunction, 

E(t) if> (x,t) = -if> (x,t), 



of the time-dependent dynamic invariant, 



E(t) = - 
v ; 2 



{p — 2ax — 5y 
J 2 



+ {px + ef 



(2.10) 



(2.11) 



d 



\a{t)a ] (t) +a ] (t)a(t)] , —(E) = 



with the required operator identity [36], [37], [125] : 

— + i ' 1 [E,H] = 0, 



1 



H = -(p 2 + x 2 ) 



The time-dependent annihilation a (t) and creation a? (t) operators are given by 



a(t)-- 
aHt) 



1 

i 

: 71 



fix + e + i 



(3x + e — i 



P 



2ax — 5 



V 



P 

2ax 



(2.12) 



(2.13) 



with p 



:-l 



d/dx in terms of solutions fl2.4p - fl2.9l) [95] . These operators satisfy the canonical 



commutation relation, 

a(t)a j (t) -a j (t)a(t) = 1, (2.14) 

and the oscillator-type spectrum of the dynamic invariant E can be obtained in a standard way by 
using the Heisenberg-Weyl algebra of the rasing and lowering operators (a "second quantization" , 
the Fock states [95J ) . In particular, 

a{t)^ {x,t) = 0, ip (x,t) = e h(t) ^ (x, t) , (2.15) 
with <p (t) = —7 (t) being the nontrivial Lewis phase [92], [125] . 
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This form of quadratic dynamic invariant and the corresponding creation and annihilation oper- 
ators for the generalized harmonic oscillators have been introduced recently in Ref. [125] (see also 
[27] . |143] and the references therein for important special cases). An application to the electro- 
magnetic-field quantization and a generalization of the coherent states are discussed in Refs. [82] 
(see also section 8) and [88], respectively. 

The key ingredients, the maximum kinematical invariance groups of the free particle and harmonic 
oscillator, were introduced in [3], @], [57], [70], [TT3], and pH] (see also [13, [76], [109], [122] . 
[141] . [142] and the references therein). We use the connection with the Ermakov-type system, 
which allows us to bypass a complexity of the traditional Lie algebra approach [91], [95] (see 
[4"5] . [8"§] and the references therein regarding the Ermakov equation). (A general procedure of 
obtaining new solutions by acting on any set of given ones by enveloping algebra of generators 
of the Heisenberg-Weyl group is described in [37]; see also [5], [7J, [3JJ], and [105] regarding the 
corresponding wavefunctions.) Finally, it is worth noting that the maximal invariance group of the 
generalized driven harmonic oscillators is isomorphic to the Schrodinger group of the free particle 
[91], [95], [113], and [TT4] . 



3. The Uncertainty Relation and Squeezing 

A quantum state is said to be "squeezed" if its oscillating variances ((Ap) 2 ) and ((Ax) 2 ) become 
smaller than the variances of the "static" vacuum state ({ApY) = ((AxY) = 1/2 (with H = 1). For 
the harmonic oscillator, the product of the variances attains a minimum value only at the instances 
when one variance is a minimum and the other is a maximum. If the minimum value of the product 
is equal to 1/4, then the state is called a minimum-uncertainty squeezed state (see, for example, 
@D, [EE], [133], [138], [139], [153], and [158]). This property can be easily verified for solution (Q. 

According to (I2.13p . the corresponding expectation values oscillate sinusoidally in time 

1 d 
(x) = -—[(2a e - p 6 ) sin t + e cost] , — (x) = (p), (3.1) 

(p) = -^-[(2a e - /Vo) cost- e sin t] , — (p) = -(x) (3.2) 

with the initial data (x)\ t=0 = —eo/(3 and {p)\ t=0 = — (2«o£o ~ Po$o) //V This provides a classical 
interpretation of the "hidden" parameters. 

The expectation values (x) and (p) satisfy the classical equation for harmonic motion, y" + y = 0, 
with the total "classical mechanical energy" given by 



(3.3) 

t=o 



For the standard deviations on solution (12.2p - (12.9l) . one gets 

((ApY) = (p 2 ) - (p) 2 (3.4) 
1 + Aal + /3q + (4a 2 , + - l) cos 2t - 4a sin It 

((Ax) 2 ) = (x 2 ) - (x) 2 (3.5) 



and 
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1 + Aa\ + /9q - (4«o + Po ~ 1) cos 2t + 4a sin 2t 



((Ap) 2 )((Ax) 2 ) 



16/? 4 



(3.6) 



x 



Here, 



(l + 4a 2 , + /3 4 ) 2 - ((4a 2 , + /?o - l) cos2t - 4a sin2t) 1 

4a 2 + /3 4 



((A P r> 



2/3 z 



a* = ((Ax)') 



2/3 



2 » 



(3.7) 



px 



^(ApAx + AxAp) = 



with two invariants: 



cr„ + <7 X . 



4a 2 + /3 4 + 1 4a 2 + 0* + 1 
■ " 2/3 2 ' 



cr 



1 

4 



2/3'' 2/31 °p* a - * " VX 4 (3 ' 8) 

(More invariants are given by in (I6.lip - fl6.12p .) A family of minimum-uncertainty states, when 
((Ap) 2 ) = ({Ax) 2 ) = 1/2, is defined by taking a = and /3 2 = 1. 

By adding fl3.3p - fl3.5l) . we arrive at 



(H) = \ [(P 2 ) + (x 2 )} 

_ 1 + Aa\ + /3 4 (2a g - /3 £ ) 2 + g 2 1 
+ 2/3 2 -2 



(3.9) 



4/3 ^ 

for the total "quantum mechanical energy" in terms of our "hidden" parameters/integrals of motion 
(the vacuum value 1/2 occurs when f3 = ±1 and «o = <5o = £o = 0). See also [15] and [59] . 

Therefore, the upper and lower bound in the Heisenberg uncertainty relation are given by 



max [((Ap) 2 } ((Ax) 



4 x2 



(l + 4a 2 + /3 4 ) 



16/3 



when cot 2t 



Aa 



and 



mm 



[((Ap) 2 )((Ax) 2 }} 



1 

4 : 



if tan 2t 



Aal + /3 4 
4a 



Aa 2 + /3o 



(3.10) 



(3.11) 



Our explicit formulas f l3.4p -f l33|) show that the product of the variances attains the minimum value 
1/4 only at the instances that one variance is a minimum and the other is a maximum as stated 
in [61]. The corresponding squeezing of one of the variances is also explicitly described by our 
formulas. Indeed, one gets 

1/2 / oX 1/2 

~ ! ^ , (3.12) 



(4a 2 + /3 4 -l)cos2t-4a sin2t = ±(^4a 2 + (/3 2 + l) 2 ) (4a 2 + (/3 2 - l) : 
under the minimization condition (13. lip and at the minimum 

1 /2 

1 + 4a 2 + /3 4 ± (4a 2 + (/3 2 + l) 2 ) (4a 2 + (/3 2 - l) 



((Ap) 2 ) 



4Z3 2 , 



1/2' 



(3.13) 
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1 + 4og + ft T (4c^ + (ft + I) 2 )' 72 (4^ + (ft - if 



(3.14) 



4/^ 

for all real values of our parameters. At this instant the squeezing occur: 

<(AP) 2 > > \ (< i) , <(A*f > < i (> i 

(for upper and lower signs, respectively). As a result, we present the minimum- uncertainty squeezed 
states for the simple harmonic oscillator in closed form fl2.4p -f l2l3]) (see also [61] for numerical 
simulations). These states form a six-parameter family (a natural generalization will be discussed 
in the next section). The corresponding wave function in the momentum representation is derived 
by the (inverse) Fourier transform of f)2.2p and fl2.3p - fl2.9l) (95]. Experimentally observed time- 
oscillations of the velocity variance |11U] reveal certain damping, which can be explain in models of 
quantum damped oscillators discussed in Refs. [25], [27], and [SB] (see also the references therein). 

Example. In a special case, one simplifies 

<(Ap) 2 ) = (P 2 ) - (V? = 1 " 2 °°, 5in2f , (3.15) 

( (A^)^VW 2 = 1 + 2 £ Sm2 ' . (3-16) 

provided that 4a 2 , + f3 = 1. In the case of the Schrodinger ground state "static" solution [131j . 
when a = 5 = e = and (3 = ±1, we arrive at (x) = (p) = and 

((Ap) 2 ) = ((Ax) 2 ) = \ (3.17) 

as presented in the textbooks [SI], [53], [3B], [HI], [HS], [108] . (In general, dependence on the 
quantum number n, which disappears from the Ehrenfest theorem [33], |60j . is coming back at the 
level of the higher moments of distribution [95].) 

According to f l3~T5|) - fl3TT6|) . 

((Ap) 2 )((Ax) 2 ) = 1 " 4 "° 4 Sin22t , 4« 2 + ^ = l (3.18) 

and the product is equal to 1/4, if sin 2 2t = 1. (For the coherent states «o = and /3 = ±1, which 
describes a two-parameter family with the initial data (x)\ t=0 = =F^o and (p)\ t=0 = ±^o.) 

Our formulas (I3.15p -f l3.16p show that once again the product of the variances attains the minimum 
value 1/4 only at the instances when one variance is a minimum and the other is a maximum [61J, 
|158j . The corresponding squeezing of one of the variances is also explicitly described. For example, 
if sin2t = 1, 

it\ \2\ l-2« 1 if \ \2\ l + 2a 1 

<(Ap) > = ~W 5' <( ' > = ~W 2 ( ' 

provided that < «o < 1/2 and 4a 2 , + /3 = 1. 
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4. An Extension: the TCS States 



We construct an analog of the coherent states (generalized coherent, or the TCS states in the 
terminology of Ref. |158j ) in a standard fashion 



oo 

i2 , 



^(x,t) = e-^/ 2 J2^n(x,t) ±= (4.1) 
r, ynl 

n=0 v 



n=0 

where ( is an arbitrary complex parameter and the "dynamic" wave functions are given by equations 



Tl I 1 V 

ex ] 

[1.2) and (1.16) of Ref. [95]: 



i(«z 2 +<5:E+ft)+i(2n+l)7 

^ n (x,t) = ^_ e-« 2 / 2 tf n (0, £ = /3x + e (4.2) 

(see also [36] and |105j ). where H n (x) are the Hermite polynomials |115j . In the explicit form [131] . 

$ U t ) = ^^ e -(€ 2 +M 2 )/2 e *K+fe+ K +7) Y(*Y ^1 (4.3) 



n=0 



1 p (r ? 2 -|r?| 2 )/2 p i(q a; 2 +<5 3; +/c+7) ff -(g-v^> 7 ) 2 /2 

and the eigenvalue problem is given by [158] : 

a (t) ip (x, t) = r]ip (x,t) . (4.4) 
An elementary calculation shows that on these "dynamic coherent states" , 

(x) =J-^( T 1 + V*)- |, <*> \t=o = y Q ICI cos (2 ( 7o + 0)) - (4.5) 



and 



/5 / *s ay/2 . /. 2ae\ 

(p) = (77 - if) + (I/ + /f) + ^ - — J , (4.6) 



(P>l 



t=o 



A, V2 |C| sin (2 ( 7o + 0)) + 2 3 / 2 ^ |C| cos (2 ( 7o + 0)) + ^ - ^) 



if C = ICI e21 ^- Moreover, a direct Mathematica verification shows that these expectation values 
satisfy the required classical equation for simple harmonic motion. 

A similar calculation reveals that the corresponding oscillating variances ({ApY) and ({AxY) 
coincide with those for the "dynamic vacuum states" given by (I3.4p - fl3.5p . The "dynamic coherent 
states" (14. 3 p are also the minimum-uncertainty squeezed states but they are not eigenfunctions of 
the time-dependent dynamic invariant (12. lip when rj ^ 0. 



The Wigner function [65], [91], [[29], [155] 

1 

2^ 



W {x, p) = / V* (x + y/2) ip{x- y/2) e vpy dy, (4.7) 



MINIMUM-UNCERTAINTY SQUEEZED STATES 



9 



for the TCS states (14. 3[) is given by 

1 



W (x,p) 



where 



P 



-[P + i 



exp 



p — 2ax — 5 



7] — T) 

~7T 



Q - 



77 + 77* 

~7T 



ft 



Q = f3x + e. 



In view of (I4.5I) - (I4.6I) . we arrive at the following expression of the Wigner function: 



W (x, p) = — exp 



7T 



(p- (p)Y 

(3 2 



Aa Aa 2 + /3 4 

+ ^2 (p - (p>) - — n2 — ( x - ( x )y 



(4.8) 
(4.9) 

(4.10) 



in terms of the classical trajectories (x) and (p) and solutions of the Ermakov-type system (12.4]) - 
(12. 5p provided /x /9 — 1- Taking into account the time-dependent variances (13.71) . one gets [31], [59] . 
□32]: 

W fop) = ~ exp [-2 (p - (p)) 2 - 2^ (p - (p» (x - (x)) + <x p (x - (x» 2 )] , (4.11) 

where a p , a x , and a px are given by ( 13. 71) . Then 

(x,p; t) = W (x cost — psmt, x sin t + pcost; t — 0) (4-12) 

by a direct calculation — the graph of Wigner function rotates in the phase plane without changing 
its shape [137] . (In a traditional approach, the quantum Liouville equation of motion for Wigner 
function of the corresponding quadratic system is used in order to determine this time evolution 
[129J. We have obtained the same result directly from the wave functions.) Some Mathematica 
animations can be found in Ref. l83l. 



5. The Moyal Functions 
The total energy of a "dynamic harmonic state" (\A.2\i can be presented as 



(H) = 2 [(P 2 ) + (* 2 >] 



n 



l\ 1 + Aal + Pi (2a e - (3 5 ) 2 + e 



V 2p 2$l 
by (A.3)-(A.5) of Ref. [95]. 

The Moyal functions [111] for the "dynamic harmonic states" (14.21) : 

1 f°° 

W mn (x, p, t) = — J r m (x + y/2, t) ^ n (x - y/2, t) e™ dy 



(5.1) 



(5.2) 



can be evaluated in terms of Laguerre and Charlier polynomials in a standard way [115] . [129] : 



W mn (x, p, t) 



( -] \m 2i(n-m)j 

[Zll e „-Q 2 -P 2 /P 2 r,(m-n)/2 



71 




(5.3) 



x Q 



iP 



2 Q 2 + 



13 ) ~ m KV'fi, 

Once again, the time evolution of the corresponding Wigner function W nn (x,p,t) is defined by 
equation (14.121) . 
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In the case of an arbitrary linear combination, 

V> (x, t) = ^ c m ip m (x, t) , (5.4) 

the Wigner function can be obtain as a double sum of the Moyal functions: 

(x,p,t). (5.5) 

771,71 

A coherent superposition of two states with n = and n = 1 was experimentally realized in 
Ref. [EED]. Moreover, the state of the electromagnetic field can be chosen anywhere between the 
single-photon and squeezed state in [72]. 



6. ElGENFUNCTION EXPANSIONS 



Experimentally observed statistics for various squeezed states of photons and ions in a box [15J, 
[59] . [90] . [102] . [106] , [128] can be naturally explained in terms of explicit developments with respect 
to the Fock states. For a linear harmonic oscillator in the coordinate representation, we discuss a 
general case by using the wave functions and known expansions in Hermite polynomials [85], [93J, 
[115] . Group-theoretical properties are discussed elsewhere. 



6.1. Familiar Expansions. For the stationary harmonic oscillator wave functions, 

-x 2 /2 

* n {x) = -?==H n (x), (6.1) 
there are two well known expansions: 

oo 

e^ r+Bx ^ n ( X + A) = J2 T mn (A B, T) m m (x) , (6.2) 

777=0 

where 

/oo 
m* m {x)e^ +Bx ^ n {x + A) dx (6.3) 



-oo 

•777 — 71 



e i(T-AB/2) g-f/2 ( ^_+_^^ ~ B 



y/m\n\ \ V% J \ V2 

x 2-F0 I -n, -m; -- 1 



v 



with v = (A 2 + B 2 ) /2 (see, for example, [93], [115] for relations with the Heisenberg-Weyl group, 
Charlier polynomials, and Poisson distribution) and 

oo 

e iax2 V n (0x) = M mn (a, 0) # m (x) . (6.4) 

777 =0 

By the orthogonality, 

POO 

M mn (a, 13)= m* m (x) e iax2 V n (fix) dx, (6.5) 
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and one can use the integral evaluated by Bailey: 



A 



i A x H m (ax) H n (bx) dx 

2 m+n f m + n+ I 



(6.6) 



m+n+1 " 



( m+ ^ +1 ) (a 2 -A 2 ) m/2 (& 2 -A 2 ) n/2 




ab 



y/(a* - A 2 ) (6* - A 2 ) 

Re A 2 > 0, if m + n is even; the integral vanishes by symmetry if m + n is odd; see Refs. [6], [98] and 
the references therein for earlier works on these integrals, some of their special cases and extensions. 
As a result, 



, r , „ , 2 m+n ^ f m + n + 1 

M mn (a,(3) =i n \l^^- T 



x 



1-/3 2 



mw.-K 

m/2 



(6.7) 



+ ia 



1-/3 2 



— ia 



n/2 



i+£ 

2 

—771, 



)Q 



-n 



x 2 Fi 



(1 — m — n) 



(m+n+1) /2 



;i|l± 



2 ' " \ A /4a 2 +(/3 2 -l)' 

The terminating hypergeometric function can be transformed as follows 
—k, —n i 



'1 — k — n) 



= < 



f ( l/2) r (1/2), 
(V2) r+S 

(3/2) r (3/2) 



^1 



— r, — s 
1/2 



(6.8) 



(3/2) 



r+s 



-r, — S 
3/2 



if k = 2r, 77, = 2s, 



-C 2 , if k = 2r + l, n = 2s + l. 



It is valid in the entire complex plane; the details are given in Appendix B of [55]. The trans- 
formation ( 16. 8 p completes evaluation of the Bailey integral (I6.6P and the matrix elements ( 16. 7p in 
terms of the hypergeometric functions. (Relations with the group SU (1,1), Meixner polynomials 
|115j . and with two special cases of the negative binomial, or Pascal, distribution [85] are discussed 
elsewhere.) 



6.2. Probability Amplitudes. Expansions (16. 2p and ( 16. 4p results in 

e i(2n+l)(7-7 ) °^ 

ip n (x, t) = ^ C mn (t) ^ m (x) 

m=0 



(6.9) 
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where 



°° / 5 \ 

C mn (t) = ^2 M mk (a, 0) T kn \e, -, nj 



(6.10) 



fc=0 

00 



fc=0 

The invariants are 

..2 «4 , 1 a„.2 , #4 



Jmfc (-0,0 — , k - — ) M kn (a, (3) . 

h — n \P 



(3 ' /3 2 



4a 2 + /T + 1 _ Aa 2 + gg + 1 _ 5e_ 5 e 

2 5 2 2 5l ; 2 / 2ae\ 2 e 2 f 2a e ^ 2 



J 2 °"/3 2 ' /3 2 
by the direct calculation. Another useful identity is given by 

4a 2 + /3 4 + l 4a 2 +(/3 2 ±l) 2 

o ± 1 = o = <7 D + (X,. ± 1. 

2/3 2 2/3 2 p 

Thus all arguments of the hypergeometric functions in ( 16.101) are constants. Moreover, the 
dependencies of the matrix elements are given only by complex phase factors: 

T mn (e, |, ft) = e 2i(m-n)( 7 - 7o ) ^ ^ ^ ^ 

/ e 2a£ ae 2 



J(n-m)t rp f £ r 2a g a 6 Q \ 



f P 0' 
and 

M mn (a,P) = e -^ m+l ^-^e-^ n+l ^\[^M mn (a ,/3 ) 

V ^0 

in view of the following identities 

6 ^ . (60 

j +lE = U 



+ is = I + ie j e 2 ^"^, 



2ae e / 2a e . .£0 , 



0/3 V A> 00 



-it 



1 — /3 2 /l — /3 2 \ 

+ ia = e~ lt [ — - — - + ia J / (2a sin t + cost + sin t) 

— iao^ / (2a sin t + cos t + sin t) 



2 

?2 



6.11) 
6.12) 

6.13) 
time- 

6.14) 
6.15) 

6.16) 

6.17) 
6.18) 
6.19) 
6.20) 



and some of their complex conjugates (see also Appendix A for a complex parametrization of the 
Schrodinger group). 

Finally, the eigenfunction expansion takes the form 

00 

V>„ {x, t) = -^z Yl °mn e-^W V m (x) , (6.21) 



^0 m=0 
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where the time-independent coefficients are explicitly given by 

Cmn = (otQ, O ) Tkn ( £q, 7p, ^0 J 



(6.22) 



oo 




2a £ 




) 



El 



mk 



Kq — 



M kn (a ,/3 ) 



k=0 



in terms of the initial data/integrals of motion (of the corresponding Ermakov-type system). The 
total probability amplitude is a product of two infinite matrices related to the Poisson and Pascal 
distributions. 

Moreover, a combination of (14. ip and (16.211) gives the eigenfunction expansion of the TCS states. 
It is worth noting also that our expansion (16.211) gives an independent verification of the fact that 
the "missing" solutions (14 .2p do satisfy the time-dependent Schrodinger equation (12.11) . Indeed, 
they are written as an explicit superposition of the standard solutions. 

7. NONCLASSICAL HARMONIC STATES OF MOTION AND PHOTON STATISTICS 

A fundamental manifestation of the interaction between an atom and a field mode at resonance 
in an ideal cavity is the Rabi oscillations [51]. The first observation of the nonclassical radiation 
field of a micromaser is reported in p_20j (the statistical and discrete nature of the photon field 
leads to collapse and revivals in the Rabi nutation |121j ). Implementation of light for purposes 
of quantum information relies on the ability to synthesize, manipulate, and characterize various 
quantum states of the electromagnetic field. A review [102] covers the latest developments in 
quantum-state tomography of optical fields and photons (see also the references therein). 



Various classes of motional states in ion traps are discussed, for example, in [90J . Our expansion 




formula (I6.22p is consistent with statistics for the coherent, squeezed, and Fock states observed in 
Refs. [15] and |106] for ions and photons in a box (see also [39] and [90]). A method to measure the 
quantum state of a harmonic oscillator through instantaneous probe-system interaction, preventing 
decoherence from disturbing the measurement, is proposed in |126j . 

7.1. Coherent States. In breakthrough experiments of the NIST group on engineering ionic states 
of motion, the coherent states of a single 9 Be + ion confined in a Paul trap were produced from the 
ground state by a spatially uniform classical driving field and by "moving standing wave" (see [90J, 
[106J and the references therein for details). For the data presented in [106J, the authors used the 
first method. The Poissonian distribution with the fitted mean quantum number n = 3.1 ±0.1 
was identified from Fourier analysis of Rabi oscillations. In our notation, ao — 0, /3 = 1, and 



Time evolution of the coherent state of cold Cs atoms was measured in |110j . For experimentally 
observed coherent photon states [5j], see, for example, [15] and [99] . 

7.2. Squeezed Vacuum and Fock States. The minimum-uncertainty squeezed state with 7 = 
= £ o — K o — is called the squeezed vacuum (see [39], [79], and [90J when a = 0). Expansion 
( 16.211) simplifies to 




n = 



^0 OM) 



(7.1) 



Vf 1 (*) v 7 ^ 
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ft,- 

+ ia 



provided that /i /?o = 1- The probability distribution is restricted to the even states and given by 

in terms of the variances (13. 8p . This is a special case of the negative binomial, or Pascal, distribution. 

A vacuum squeezed state of ionic motion was created in the NIST group experiments [106J by a 
parametric drive at 2u (see also [HI], [SD] and the references therein). The data were fitted to the 
vacuum state distribution (17. 2p with a v + a x = 40 ± 10 and «o = (corresponding to a noise level 
16 dB below the zero-point variance in the squeezed quadrature component; see [90] and [106J for 
more experimental details). 

A vacuum squeezed state of motion of neutral Cs atoms was also generated in [110] . Here, the 
cold atom sample containes about 10 5 atoms. Therefore a single image provides the full velocity 
distribution of the quantum state and the squeezing can be readily visualized — a set of images 
gives the state's time evolution (see [HOj and the references therein for more details). 

In a similar fashion, for the squeezed Fock state with n — 1 and 7 = 5$ = Eq = = 0, expansion 
(I6.2ip simplifies to 



/i (t) Vtt 



1-/3 o , • 

1 ta 



Po v 2 p+1 r( P + 3/2) V 2 V p ^2 P+m t * 

[-2- - m ° 

The corresponding Pascal distribution for the odd states is given by 

Pm=2p+1 = ( ' > (7-4) 

(a p + a x + lf /2 p\ \cT p + a x + lJ 

where (3/2) = 1 and (3/2) = (3/2) (5/2) • • • (1/2 + p) . These squeezed Fock states were generated 
in [TT] and their dynamics was analyzed in [llCJj . When e 7^ 0, displaced Fock states of the 
electromagnetic field, have been synthesized in [99] (see also the references therein). 

Moreover, even/odd oscillations in the photon number distribution of the "squeezed vacuum" 
state, which are consequence of pair-wise generation of photon, were observed in [15], |128j . For 
an ideal minimum-uncertainty squeezed state zero probabilities for odd n are expected, since the 
Hamiltonian describing the parametric process occurring inside the nonlinear crystal is quadratic 
in the creation and annihilation operators [39] ■ |129j . However, the probabilities for odd photon 
numbers are nonzero because the squeezed state detected there is a mixed state having undergone 
losses inside the resonator and during the detection process which cause the distribution to smear 
out (see [39] and |128] for more details). The corresponding Pascal distributions (17. 2p and (I7.4p 
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have different parameter values for even and odd states, which is consistent with the result of these 
experiments. Further details will be discussed elsewhere. 

7.3. Engineering Mixed Squeezed States. Generation of a coherent superposition of the ground 
state and the first excited states of motion of cold Cs atoms in the harmonic microtraps, namely, 

ip (x, t) = c ^o (x, t) + ci^i (x, t) , (7.5) 

where Co = 2~ 1//2 and Cl = 2" 1 / 2 e^, was reported in [110] and the corresponding time evolution 
had been experimentally observed. This evolution is obviously nonclassical and contrasts with that 
of a coherent state which oscillates as a classical particle without deformation (see |110j for more 
details). This dynamics is consistent with our expansion (I6.22p but details of these calculations will 
appear elsewhere. 

8. An Application to Cavity QED and Quantum Optics 

Foundations of quantum electrodynamics and quantum optics are presented in many excellent 
books and articles 0, [E], 0, [13], [2H], [33], gj, @S], @8], [IS], jSg, |M], HU, |5S], US], 
PI, PI, EE], H23, P23, [132], [135], [152], [153], and [151]. Here, we suggest a modification of 
the radiation field operators in a perfect cavity in order to incorporate the Schrodinger symmetry 
group into the second quantization. Our approach gives a natural description of squeezed photons 
that can be created as a result of parametric amplification of quantum fluctuations in the dynamic 
Casimir effect [M], |156] and are registered in quantum optics [15], |102j . [116] . 

8.1. Radiation Field Quantization in a Perfect Cavity. In the formalism of second quantiza- 
tion, one expands electromagnetic fields in terms of resonant modes of the particular cavity under 
consideration [¥!] . [73] . |129j . |135j . The cavity is represented by a volume V, bounded by a closed 
surface. Let E a (r) , k\ = oj 2 a /c 2 be the eigenfunctions and the eigenvalues of the corresponding 
boundary-value problem: 

V X V X E-k 2 E = in V (8.1) 

n x E = on 5, 

where n is a unit normal vector to S. The vector functions H a (r) are related to E a (r) by 

V x E a = k a H a , V X H a = fc a £ Q . (8.2) 
The eigenfunctions are orthonormal in V : 

/ E a -E p dV = 5 a p, [ H a -H p dV = <V (8.3) 

The electric and magnetic fields are expanded in the following forms 

E (r, t) = -Va^J^P- (*) ^« ( r ) ' ( 8 - 4 ) 

a 

(t) H a (r) . 

a 

The total energy is given by 

H= r dv= ^{pl + colql) (8.5) 
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and the Maxwell equations, 

vx.af, vxir = if, (8.6) 

c at c at 

are equivalent to the canonical Hamiltonian equations, 

— - — -p ^--^H-- L0 2 q (8 7) 

dt dp a a ' dt dq a a 

respectively. 

In the second quantization, one replaces canonically conjugate coordinates and momenta by 
time-dependent operators q a (t) and p a (t) that satisfy the commutation rules 

[q a (t) , qf) (*)] = [p a (t) , p p {t)\ = 0, [q a (t) , Pfj {t)\ = ih5 a p. (8.8) 

The time-evolution is determined by the Heisenberg equations of motion |60j : 



j t Pa (t) = ^ [ Pa (t) , H) , j f q a (t) = % - [q a (t) , H) , (8.9) 

with appropriate initial conditions^ (From now on, we consider a single photon cavity mode, say 
a, with frequency co a = 1 and use the units c = h = 1.) 

8.2. Nonstandard Solutions of Heisenberg's Equations. Explicit solution of equations (18. 9|) 
for squeezed states can be found as follows 

p(t) = !fiL|M , W = !M+M. ( 8 ,o, 

The time-dependent annihilation b (t) and creation (t) operators are given by [S2] 

b(t) = — \^x + e + i P - J , (8.11) 

r+ / \ e 2i7 / p - 2ax - 5 

tf(t) = —=[Px + e-i F 



V2\ P 

in terms of solutions (I2.4p -( l2~9|) of the corresponding Ermakov-type system. The time-independent 
operators x and p obey the canonical commutation rule [x,p] = i in an abstract Hilbert space. At 
all times, 

6 (*)#(*) -6 f (t)6(t) = 1. (8.12) 
By back substitution, operators b (t) and (t) are solutions of the Heisenberg equation: 

d- 
dt' 

with the standard Hamiltonian 



-b(t) =i\b(t), H], H ] (t) = i\ ? (t) , h] , (8.13) 



= I (p 2 + x 2 ) (8.14) 



subject to the following initial conditions 



b (0) = —j=- I (3 x + e + z -g- ), (8.15) 



2 The standard form of Heisenberg's equations can be obtained by the time reversal t — > — t (with ao — > —&0i 
7 — > — 7 , 5q — > —So, and kq — > —kq; see below). 
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&t( ) 



e 2 ^o ( p _2a x-5o 
p x + e — i— 



V2 



o 



The creation and annihilation operators (18. lip allow us to incorporate the Schrodinger group of 
harmonic oscillator, originally found in the coordinate representation [114] . into a more abstract 
Heisenberg picture — the classical case occurs when (3 = 1 and a = 7 = 5 = e = n = 0. (For 
the sake of simplicity, we have restricted ourself to a single photon mode with frequency uj = 1.) 



8.3. Dynamic Fock Space for a Single Mode. The time-dependent quadratic invariant, 

[p — 2ax — 5) 2 



E(t) 



with 



dE 



2 



+ (fix + sY 



8.16) 



d 
dt 



(E{t)) = 



E,H 



1 



0, H = -(p 2 + x 2 ) 



(8.17) 



extends the standard Hamiltonian/Number operator H for any given real values of the "hidden" 
parameters/integrals of motion in our description of the squeezed photon state. The oscillator- type 
spectrum, 

r. , 

~.i< 



(8.19) 



E{t)\il> n (t))=\n+-J |V n (f)>, 
can be obtained by using the modified creation and annihilation operators [2]: 

6(*)K(*)) = V^ |^n-l(*)>, 

?(t)i^ n (t)) = v / ^n \^ n+1 (t)). 

For the "minimum-uncertainty squeezed states" , one gets 

b(t)\% (t)) = 

with 



(Vo (t)\H\%(t)) 



+ 



> 



4/9^ 2/3^ " 2 

in the Schrodinger picture. The generalized coherent (or TCS's) states are given by 

&(*) |V> (f)> = CliK*)> 

for an arbitrary complex ( ^ 0. 



(8.20) 
(8.21) 

(8.22) 



8.4. Expectation Values and Variances for Field Oscillators. The noncommuting electric 
E(r,t) and magnetic H (r,t) field operators are given by equations (18.41) and (I8.10p - fl8.lip for a 
squeezed photon in the Heisenberg picture, which provides a more direct analogy between quantum 
and classical physics [59]. The electromagnetic radiation mode in a cavity resonator is analogous 
to a harmonic oscillator [61]. In the Schrodinger picture, all previous results on the minimum- 
uncertainty squeezed states can be reproduced for the field oscillators in an operator QED-style. 
For a single mode with u a = 1, 

(E (r, t)> = - V^rX (r) (*) \p\ ^ (*)> , (8.23) 
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// (r, t)> = V^H a (r) (ij n (t) \x\ ^ n (t)> , 

where equations (I3.ip - fl3.2p hold. The corresponding variances are given (up to a normalization) 
by equations (A.4)-(A.5) of Ref. [95]. 

The minimum-uncertainty squeezed states are identified in quantum optics [30j , [61], [56], [72] 
[90], [133], [136], [IIS], [ES], USE] and in state tomography [T6], [H], [SI], [102]. They are also 
important in the dynamical Casimir effect [31], [32], [33], [H], [52], [82], [M], [M], [156], and [T57J, 
where the photon squeezing occurs as a result of a "parametric excitation" of vacuum oscillations. 



9. An Important Variance 



The Hamiltonian H = (p 2 + x 2 ) /2 can be rewritten in terms of the creation and annihilation 
operators (I2.13P as follows: 



H 



Aa 2 - (3 4 + 1 
A~f 



— ia I a 2 (t) + 



Aa 2 - /3 4 + 1 

4/f 



+ ia ] a) (t) ' 



(9.1) 



Aa 2 + (3 4 + 1 r ^ , . _ + . . ^ , . ^ . . -i 
+ [a (t) a j (t) + a ] (t) a (*)] 



+V/2 
+V2 



4/3 








+ 2 U 



2ae 

T 



2ae 

T 

2ae 

+ 



e 


i/3 


(- 


2ae 


2/3 2 


~ ~2 





e 
2(3 2 


i/3 

+ T 


(- 


2ae 





a(t) 
a ] (t) 



2/3' 



and by definition: 

Var H = ((H - (H)f) = (H 2 ) - (H) 
Then a direct Mathematica calculation results in 



Var H 



Aa 2 +% + \) 2 ) (4ag + - 1) 



n 



+ 



(Aa 2 + (3t + l) ((2a e -(3 5o) 2 + e 2 ) 







2 , 4 

2 








■ 3 

+ i 

n + - 



(9.2) 



(9.3) 



for the wave functions (14. 2p in terms of the invariants f l6.lip -f l6.12l) . (These calculations can be 
performed in pure operator form with the help of standard relations (1.15) of Ref. [95]; see also 
f)8.19p .) In terms of the variances, 

2 



Vartf = l - [K + a x f - l] 



1 

n H — 
2 



(9.4) 



+2 [a p (p) 2 + 2a px (p)(x) + a x (x) 2 ] (n+~ 



where a p , a x , and a px are given by f l3.7p . When n = 0, this formula is consistent with the variance 
of the number operator derived for a generic Gaussian Wigner function in Ref. [39J. 
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A similar expression holds for the TCS states. Computational details are left to the reader. 



10. A Conclusion 



In this paper, we review the nonclassical states of harmonic motion, which were originally found 
in [1U5] (in the coordinate representation) and have been rediscovered recently in m a form 
which is convenient in applications to cavity QED and quantum optics. In particular, the minimum- 
uncertainty squeezed states are studied in detail. Expansions in the Fock states are established and 
their relations with experimentally observed photon statistics are briefly discussed. In the method of 
second quantization, a modification of the radiation field operators for squeezed photons in a perfect 
cavity is suggested with the help of a nonstandard solution of Heisenberg's equation of motion. 
These results should be of interest to everyone who studies introductory quantum mechanics and 
quantum optics. 
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Appendix A. A Complex Parametrization of the Schrodinger Group 



Let us introduce the following complex-valued function: 

z = Cl e il + c 2 e~ i \ z" + z = 0, (A.l) 

where 

ci = (1 + Pi) /2 - ia , c 2 = (l-#j)/2 + za (A.2) 
(ci + c 2 = l, |ci| 2 - |c 2 | 2 = Pi) , 

and 

c 3 = - ie . (A. 3) 

Po 

Then equations ( 12 . 3 H — ( T2T9T) can be rewritten in a compact form in terms of our complex parameters 
ci, c 2 , and C3. With the help of identities (I6.17l) -( l6.20p . one gets 

^ _ 11 _ (\„ |2 , „ „*„2it , „*„ „-2it , |„ |2\ 1/2 
t*0 



| Cl | 2 + Cl cle 2U + c\c 2 e~ m + |c 2 | 2 ) 1 (A.4) 



and 



.cic* 2 e 2it - c\c 2 e 2it . r , 

01 = 1 2VT 2 ' 



B n \l \c\\ 2 - \c 2 \ 2 

P = pr = ± V , , , (A.6) 

\z\ \z\ 

7 = 7 -iargz, (A.7) 
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— (c 3 e iaigz + c* 3 e~ iaisz ) , (A.l 



2\z 
i 



-(c 3 e^ z -c* 3 e~^ z ), (A.9) 

« = «0"iK(l- - c 3 2 (1 - e" 2iar ^)] . (A.10) 

o 

The inverse relations between the essential, real and complex, parameters are given by 



«o = ^ {c\c* 2 - c*c 2 ) , p = ±y|ci| 2 - |c 2 | 2 , (A.11) 

#o = ±^|ci| 2 - |c 2 | 2 (c 3 + 4) , e = l(c 3 -c;). (A.12) 

These formulas (1A.5j) - (1A.10j) provide a complex parametrization of the Schrodinger group of "hid- 
den" symmetry for the simple harmonic oscillator found in Ref. |114j (see also the explicit action of 
this group in [95], specifically equation (32) there). A similar parametrization for the wave functions 
( 14. 2 p was used in Ref. [35] (see also 
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